VASP输入文件之INCAR

您所在的位置:网站首页 vasp neb计算 VASP输入文件之INCAR

VASP输入文件之INCAR

2023-03-16 03:20| 来源: 网络整理| 查看: 265

1、输入文件介绍

简介:VASP程序的主要输入文件有INCAR,POTCAR,POSCAR和KPOINTS。INCAR文件控制了vasp进行何种性质的计算,POTCAR文件包含了体系中各类元素的赝势,POSCAR文件描述了所计算的体系的晶胞参数(包括基矢或平移矢量,晶格常数,原子位置等信息),KPOINTS描述了不可约布里渊区中k点取样,即k点设置。

2、INCAR文件关键词2.1、关键词分类对所计算的体系进行注释:SYSTEM;定义如何输入或者构造初始的电荷密度和波函数:ISTARTICHARG,INIWAV;定义电子的优化:①平面波切断动能和缀加电荷时的切断值: ENCUT,ENAUG;②电子部分优化的方法: ALGO,IALGO,LDIAG;③电荷密度混合的方法:IMIX,AMIX,AMIN,BMIX,AMIXMAG,BMIXMAG,WCNIMIX,MIXPRE,MAXMIX;④自洽迭代步数和收敛标准: NELM,NELMIN,NELMDL,EDIPE。定义离子或原子的优化:①原子位置优化的方法、移动的步长和步数:IBRION,NFREE,POTIM,NSW;②分子动力学相关参数:SMASS,TEBEG,TEEND,POMASS,NBLOCK,KBLOCK,PSTRESS;③离子驰豫收敛标准:EDIFFG。定义态密度积分的方法和参数:①smearing方法和参数;ISMEARSIGMA;②计算态密度时能量范围和点数:EMIN,EMAX,NEDOS;③计算分波态密度的参数: RWIGS,LORBIT。其他:①计算精度控制:PREC;②磁性计算:ISPIN,MAGMOM,NUPDOWN;③交换关联函数:GGA.VOSKOWN;③计算ELF和总的局域势:LELF,LVTOT;④结构优化参数:ISIE等等。其他为列出的,请参见VASP的英文manual。VASP是如何判断收敛的?可以参考:竹一:VASP是如何判断收敛的?2.2、常用关键字详解SYSYTEM:给计算体系起名字,可以随便写,起一个注释的作用。ISTART:可取0或1或2或3。①ISTART = 0意思根据INIWAV初始化波函数;②ISTART = 1意思是从WAVECAR(用于储存波函数信息)读取波函数(重定义平面波);如果没有WAVECAR文件或者WAVECAR里信息和当前计算不兼容,则波函数由程序初猜产生( ISTART= 0 ); ③ISTART = 2只在计算体积-能量相关曲线的时候用到;④ISTART = 3在重启分子动力学计算的时候使用,同时需要WAVECAR、CHGCAR(用于储存电荷密度信息)文件。建议:用ISTART = 1 是最稳妥的方法,除非特殊说明!ICHARG:可取0或1或2或4或10或11或12。①ICHARG = 0代表由WAVECAR计算电荷密度信息,若没有WAVECAR文件,则默认ICHARG = 2;②ICHARG = 1代表由CHGCAR文件读取。如果没有CHGCAR文件,则由初猜产生③ICHARG = 2代表由程序初猜产生④ICHARG = 4在重启分子动力学计算时使用⑤ICHARG = 11/12/13在做非自洽计算时使用、DOS和能带计算时使用。可参考:ICHARG - Vaspwiki。建议:用ICHARG = 1 是最稳妥的,除非说明做DOS和band struct。ENCUT:平面波截断能。ENCUT = largest ENMAX on the POTCAR file;解释:由三个正弦函数叠加可以粗糙地得到想要的波函数。要想更精确的描述波函数,显然平面波的基组越大越好,但同时计算量也会增加。ENCUT就是控制使用平面波基组大小的参数,即平面波截断能。ISMEAR:可取-5 / -4 / -3 / -2 / -1 / 0 / [integer] > 0;默认值为1;建议:金属 ISMEAR = 1 ;SIGMA = 0.2;绝缘体或半导体 ISMEAR = 0 ;SIGMA = 0.05;分子计算: ISMEAR = 0 ;SIGMA = 0.01 ISMEAR - Vaspwiki。PREC:计算精度控制;PREC = Low | Medium | High | Normal | Single | Accurate ;详情可参考:PREC - VaspwikiEDIFF和EDIFFG:收敛判据。参考2.1。3、实例讲解

下面看Si的INCAR中涉及的相关参数:

General: SYSTEM = fcc Si !自洽计算fcc结构的Si。 ISTART = 0 !开始新的计算。 ICHARG = 2 !从原子的电荷密度重叠构造初始电荷密度。 ENCUT = 240 !截断能。 ISMEAR = 0; SIGMA = 0.1 !采用Gaussian smearing方法,展宽为0.1eV。 PREC = Accurate !计算精度。



【本文地址】


今日新闻


推荐新闻


CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3